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ABSTRACT 

We present results of our search for X-ray line emission associated with the radiative decay of 
the sterile neutrino, a well-motivated dark matter candidate, in Suzaku Observatory spectra of the 
Ursa Minor dwarf spheroidal galaxy. These data represent the first deep observation of one of these 
extreme mass-to-light systems and the first dedicated dark matter search using an X-ray telescope. 
No such emission line is positively detected, and we place new constraints on the combination of the 
sterile neutrino mass, mgt, and the active-sterile neutrino oscillation mixing angle, 6. Line flux upper 
limits are derived using a maximum-likelihood-based approach that, along with the lack of intrinsic 
X-ray emission, enables us to minimize systematics and account for those that remain. The limits 
we derive match or approach the best previous results over the entire 1-20 keV mass range from a 
single Suzaku observation. These are used to place constraints on the existence of sterile neutrinos 
with given parameters in the general case and in the case where they are assumed to constitute all of 
the dark matter. The range allowed implies that sterile neutrinos remain a viable candidate to make 
up some - or all - of the dark matter and also explain pulsar kicks and various other astrophysical 
phenomena. 

Subject headings: dark matter - galaxies: dwarf - galaxies: individual (Ursa Minor) 



1. INTRODUCTION 

The observational evidence that nonbaryonic dark 
matter comprises most of the mass in the universe is 
extremely strong, but the nature of dark-matter particles 
remains a mystery. No particle included in the Standard 
Model of particle physics has the characteristics required 
to explain the dark matter. In its original formulation, 
the Standard Model described all known particles, 
including neutrinos, which were assumed to be massless. 
The discovery that neutrinos have mass has forced one 
to go beyond this minimal model: a modest extension 
of the Standard Model by gauge singlet fermions can 
accommodate the neutrino masses. One or more of 
these gauge-singlet fermions can have Majorana masses 
below the electroweak scale, in which case they appear 
as sterile neutrinos in the low-energy theory. If one of 
the sterile neutrinos has mass in the 1-20 keV range 
and has small mixing angles with the active neutrinos 
(as expected), s uch a particle is a plausib le candidate 
for dark matter (jDodelson fc Widrowlll99^ . The same 
particle could be produced in a supernova explosion, and 
its emission from a cooling neutron star could explain 
the pulsar kicks (K uscnko & Seerc 1997; Fuller ct alj 
120031 ; iKusenkol [20041 : IKusenko. Mandal. fc Mukheriei 
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l2008f l. could facilitat e core collap se su- 
pernova exp l osions (jFrver fc KusenkoT 120061 ; 
IHidaka fc Fulieil [200l and can affect the forma- 
tion of the first stars fBiermann fc Kusenkol [20061 ; 
IStasielak. Biermann. fc K uscnko 2007)_ and black 



holes (jMunvaneza fc BiermannI 120051 . 1200l . Therefore, 



there is a strong motivation to search for signatures of 
sterile neutrinos in this mass range. 

The most promising way to discover (or rule out) 
relic sterile neutrinos is with the use of X-ray telescopes 
(XRTs), su ch as Suzaku. The sterile neutrinos can deca; 



radiatively (jPal fc Wolfensteinlll98^ iBarger et ai][l9& 
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to X-ray observation (lAbazaiian. 


Ful 


er. & Patel' '20011: 


Abazaiian. Fuller, fc Tuckeij 12001 




Dolffov fc Hansen 


20021). Since this is a two-bodv decav. a line is expected 



in the X-ray spectrum. The dwarf spheroidal galaxies 
are ideal targets for these observations be cause of their 
proximity, high dark matter density (jStrigari et al.l 
HoOSa), and absence of additional obscuring X-ray 
sources (see below). The Ursa Minor and Draco dwarf 
spheroidals are the optimal targets in this class for large 
field-of-view X-ray spectroscopic investigation, based on 
their average dark matter surface densities. From XMM- 
Newton observations of the Ursa Minor dwarf spheroidal, 
sterile neutrino line flux limits comparable to the best 
previous constraints were derived - even though strong 
background flaring limited the amount of good observ 
ing t ime to 7 ks (jBovarskv. Nevalainen. fc Ru chavsk: 
I2007f ). The Suzaku Observatory ( jMitsuda et al. 200 
provides the most sensitive instruments for current 
searches for weak sterile neutrino radiative decay lines 
in the ~ 0.5 — 10 keV bandpass bec ause of its low and 
stable background (|Tawa et al.ll2008[ ). and the relatively 
sharp spectral resol u tion of its CCD sp ectrometers 
(jKovama et al.l 120071 ; iNakaiima |2008[ ). In this 
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paper, we discuss our analysis of Suzaku observations 
of the Ursa Minor system, considerably improving on 
the XMM-Newton constraints; a companion article on 
the Draco system is in preparation. These represent 
the first substantial X-ray data of such extremely 
dark-matter-dominated systems. 

1.1. Context 



Sterile neutrinos may 
resonant oscillatio ns, as proposed i 
Widrow (.Dodelson fc Widrowl I19 94D. 



be produced through non- 
proposed by Dodelson and 



and 



via various 



2002; Shi & Fuller 199 



other mechanisms (iDolgovl . 

Shaposhnikov fc TkachevI l2006riKusenkQ ,2006: iKadot 
2008t iPetraki fc Kusenkol |2008[ ). A generic predic- 
tion is that relic sterile neutrinos must decay into a 
lighter neutrino and a photon. Since this is a two- 
body decay of a particle with mass mst, decay pho- 
tons produce an intrinsically narrow line with energy 
— ms t/2 that can be detected u s ing X -ray spec- 



troscopy lAbazaiian. Fuller, fc Tuckeil 1200 Ih . While 
the decay timescale, F^"^, greatly exceeds the Hub- 
ble time, the emissivity may reach detectable levels 
in regions of large dark matter concentration. Up- 
per limits on Fgt f rom X-ray observations of the 
cosmic backgr ound 
iBovarskv etall 12006 



Abazajian fc Koushiappas '2006]; 
jAbazaiian ct al. 2007), galaxy 



clust e rs (Bovarskvet aLl" 2006b: Ricmcr-Sorcnsc n et al.l 



I2007t iBovarskv. Ruc havskiv. fc Mar kevitch .2008il and 
M31 (IWatson e t al. 2006; Ricmcr-S orensen et al.l l2007t 
IBovarskv etHi 12008( 1 may be complemented and 
improved by observ ing dwarf spheroidal galaxies 
(IBovarskv et al.ll2006c[ ). 

There is no consensus as to whether dark matter is 
cold or warm based on observations of small-scale struc- 
ture. Cold and warm dark matter reproduce large-scale 
structure equally well, but differ in their predictions on 
scales comparable to the free-streaming length. Lya for- 
est clouds provide a tracer population of the pertur- 
bation spectrum on the appropriate (non-linear) scale 
that may be used to place an upper limit on the free- 
streaming length (Narayanan et al. 2000), and have been 
interpreted as implying that the mass of the dark mat- 
ter sterile neutrino, nisti must be great er than ~ 8 keV 
(|Viel et al.l 120081 : IBovarskv et"all l2009( l for production 
exclusively by the Dodelson- Widrow (DW) non-resonant 
mixing mechanism. However, A'^-body simulations and 
observ ations of dwarf spheroidal galaxies imply a lower 
mass (|Gilmore et al.l l2007t IStrigari et all HMzi)- In 
scenarios with alternative production mechanisms that 
do not involve oscillations, sterile neutrinos may effec- 
tively be colder, allowing masses as low as 3 keV to 
be consistent with the Ly a observations (Kusenkoll2006l : 
IPetraki fc Kusenkol[200alPetrakii,20081 In this case, the 
photome tric and kinematic data from dwarf spher oidal 
galaxies (|Gilmore et al.l 120071 : IStrigari et~al\ HqOT^) can 
be satisfied by dark matter in the form of s terile neutri- 
nos w ith mass between 0.5 keV and 1.3 keV (|Bovanovskvl 
|2008| ). Constraints on sterile neutrino parameters 
may be further relaxed in mixed dark m atter models 
(|Palazzo et al.ll2007t IBovarskv et"aLll2009l ) and in mod- 
els where sterile neutrinos are produced in resonant os- 
cillati ons in the presence of significant Icpton asymme- 
tries (|Shi fc Fulled Il999t ITaine fc Shaposhnikov ' ,2008) . 
Since the sterile neutrino distribution function may sub- 



stantia ll y deviate from thermal equilibr ium (|Abazaiianl 
I2006allbl : |Petrakill2008l : lBovanovskvil2008f ). simple scaling 
arguments and alterations to the power spectrum (e.g., 
imposing a cutoff) do not suffice. New large-scale struc- 
ture formation simulations with the appropriate revisions 
are required to re-evaluate the constraints from L ya data 
(|Laine fc Shaposhnikovl 120081 : lBovanovskvll2008D . More- 
over, their interpretation relies on a thorough under- 
standing of the uncertain thermal and ionization evolu- 
tion of the intergalactic medium, which may be affected 
by the propert ies and interactions of dark-matter par- 
ticlesJBicrmannJU£usenkd 120061; iGao fc Theim3 120071 : 
lMeiksin.,2009) . Hence, we present limits solely derived 
from our X-ray data. 

The sterile neutrino mass and the mixing angle 6 deter- 
mine the rate of radiative decay to active neutrinos and 
X-ray ph otons: Vs ~» The ra diative decay width is 
equal to ijPal fc Wolfensteinlll98^ 
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2567r4 ™ 
1 

1.8 X 1031s VlO-i"/ VkeVV 

1 — 5r<«'sr-2 



(1) 



where a^.^ = 1/137.0 and = 1.2 x lO-^GcV" 

For a given production mechanism, one may be able to 
eliminate one of these two parameters by requiring that 
the total density match that of dark matter in the uni- 
verse as inferred from observations. However, given the 
various possible mechanisms for producing sterile neu- 
trinos, and the possibility that sterile neutrinos m ay ac- 
count for some - but not all - dark matter ( Palazz o et al.l 
l2on7h . one should keep an open mind and search all the 
available parameter space in the mst — & plane. 

2. DATA ANALYSIS 

2.1. Observations and Initial Reduction 

Ursa Minor (OB SID = 8020 52010) was observed with 
Suzaku (|Mitsuda et all l2007f ) between 2007 April 5 
18:33:24 and 2007 April 6 8:08:14 (UT) with an on-source 
exposure time of 70.90 ks, 57.5 (13.4) ks in a 5 x 5 (3 x 3) 
editing mode (this refers to the pixel area utilized in 
event identification that is included in the telemetry; see 
Koyama et al. 2007). At the time of observation, three 
co-aligned, 17'.8x 17'.8 field-of-vi ew X-ray Imag i ng Sp ec- 
trometer (XIS) CCD cameras (|Koyama et al.l |2007( ) - 
two front-illuminated (FL XISO and XIS3) and one back- 
illuminated (BLXISl) - were operational. The BI chip is 
more sensitive below 1 keV but has a much higher inter- 
nal background, especially above 7 keV, so that BI and 
FI CCDs complement each other for broad-band studies 
such as ours. Each XIS lies in the focal plan e of an XRT 
with a 2' half-power diameter (jSerlemitsos et al.ll2007t l. 

Suzaku guest observers are provided with unfiltered 
and pre-screened event files. We initiate our data reduc- 
tion with the Ursa Minor unfiltered event files. These 
underwent Version 2.0.6.13 pipeline processing that, for 
each event, assigns a pixel quality status, calculates the 
sky coordinates using the curre nt aspect solution an d 
correcting for spacecraft wobble (|Uchivama et al.l^2008^ ■ 
and converts pulse height amplitude (PHA) to PHA In- 
variant (PI) values that map to photon energy in a one- 
to-one manner. Observations were conducted utilizing 
the space-row charge injection (SCI) technique that re- 
verses the degradation in energy resolution caused by 
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accumulated radiation damage and aff ects the conver- 
sion from the PHA channel to energy (jNakaiima et al.l 
|2008() . Version 2 pre-processing enables one to properly 
account for the eff ect of SCI on instrumen t characteristics 
and performance (jUchiv ama "eiral][2009h . We reprocess 
the unfiltered event files by hand in order to apply up- 
dated calibration data and software. Our analysis gen- 
erally follows the procedures outlined in "The Suzaku 
Data Reduction Guide," ^ as implemented in HEAsoft 
version 6.5.1^. We recalculate PI values and grades, se- 
lect event grades (0, 2, 3, 4, 6) that correspond to X- 
ray photon events, filter on pixel status (eliminating bad 
charge transfer efficiency columns, and rows invalidated 
by the charge injection process) and select good time 
intervals (GTI) based on pointing, data and telemetry 
rates, SAA proximity ("SAA_HXD = 0, T_SAA.HXD > 
436"), and proximity to the Earth's limb and illumi- 
nated Earth ("ELV > 5, DYE_ELV > 20"). In addition, 
telemetry-saturated frames and calibration source pho- 
tons are screened out and hot and flickering pixels are 
removed. Finally, we accept only GTI where the revised 
geomagnetic cut-off rigidity C0R2> 4, thus eliminat- 
ing intervals with the highest particle background level 
(|Tawa et al.ll2008l ) without compromising overall statis- 
tical accuracy (experiments with more stringent criteria 
- C0R2> 6, C0R2> 8 - yielded consistent results, but 
with significantly poorer statistics). 5x5 event files are 
converted to a 3 x 3 mode format, and merged with the 
5x5 event files. 



2.2.1. 



2.2. Spectral Analysis Procedure 
Generation of Spectra and Spectral Response 



The extraction region used to create source and par- 
ticle background spectra is constructed to include the 
entire detector fields-of-view, excluding an 8' diameter 
circular region centered on the one bright point source in 
the field. This radius was found to balance the desires 
to avoid contamination of the diffuse emission from the 
wings of the point-spread function and to maximize the 
statistical accuracy through inclusion of regions domi- 
nated by diffuse emission. 

The spectral redistribution matrices (rmf) are gener- 
ated based on the characteristics and configuration of the 
instruments at the time of observation. The rmf and 
spectral files are binned to 2048 channels. The effective 
area functions (arf) are generated via a Monte Carlo ray- 
tracing program that calculates the efficiency of detect- 
ing a photon originating from a specified source geome- 
try and surface brightness distribution wi thin a specified 
spectral extraction region (jlshisaki et al.ll2007l ). We in- 
put 2,000,000 simulation photons per energy bin from a 
uniform source of radius 20', cast the spectral extraction 
regions into detector coordinates, and adopt a single arf 
file for all components (the Ursa Minor dark matter den- 
sity is sufficiently flat over this radius that any correction 
to such a component is small). The resulting arf file is 
scaled such that the flux calculated for spectral models 
corresponds to the full 20' circular region. 

The spectra from the FI chips, XISO and XIS3, are 
co-added and a weighted XlSO-f 3 response function cal- 
culated. 

* http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
® http:/ /heasarc. gsfc.nasa.gov/docs/software/lheasoft/ 



2.2.2. Background Components 

Since dwarf spheroidal galaxies are intrinsically weak 
X-ray emitters (see below), their X-ray spectra are dom- 
inated by internal and astrophysical background compo- 
nents. The XIS background includes contributions from 
non- X-ray charged particle background (NXB), galactic 
X-ray background (GXB), and (extragalactic) cosmic X- 
ray background (CXB). 

The NXB component is estimated from observations 
of the night earth taken in the SCI mode within 150 
days of the starting or ending dates of the Ursa Minor 
observation (|Tawa et al.ll2008l i. The NXB event list in 
that time interval undergoes the identical screening as 
the source data, is sorted by geomagnetic cut-off rigid- 
ity, and weighted according to t he cut-off rigidity distri- 
bution in the source event file (|Tawa et al.l 120081 ). The 
estimated NXB spectra includes only those events col- 
lected in the region on the detector from which the source 
spectrum is extracted. XISO and XIS3 NXB spectra are 
co-added. The effective NXB spectra exposure times 
and count rates for the XIS0-h3 (XISl) are 1.1 (0.55) 
million seconds and 0.062 (0.090) counts per second, re- 
spectively, where the count rates are integrated over the 
0.45-10.4 (0.45-7.6) keV bandpass (see below). 

The GXB may be characterized by a two-temperature 
thermal plasma (corresponding to halo and local hot 
bubble contributions) and the CXB by a single power 
law. We find that this model produces a good global si- 
multaneous fit (x^ = 1205 for 1232 de grees of freedom) 
to the NXB-subtracted, binned, 0.45-10.4 keV XIS0-h3 
and 0.45-7.6 keV XISl spectra (Figures 1 and 2), where 
the two spectra are fitted by the identical model save 
for a constant multiplicative factor. Additional compo- 
nents are not required; however, allowing non-solar oxy- 
gen abundance does improve the flt. We find that the 
resulting parameters are consistent with the spectrum 
extracted from the i? 05.4 T All-Sky Survey^". 

Any intrinsic X-ray continuum is indeed negligi- 
ble. With a maximurn circular velocity < 30 km 
(jStrigari et al.l [2007aD . X-ray emitting gas cannot 
be bound to the Ursa Minor dwarf spheroidal. We 
may estimate the emission from unresolved stellar X- 
ray sources by scaling to the optical luminosity of Ursa 
Minor (~ 3 x 10^ L(7, ) using the relat i on an d X-ray 
spectral parameters in iRevnivtsev et al.l (|2008f l. There 
may also be an intermediate mass black hole in Ursa 
Minor. Based on extensions of the observed correla- 
tions of nuclear black h ole mass with velocity dispersion 
("Geb hardt et~al] 12000: Ferr arese fc FordI 120051: iGraham 
| 2008al H), bulge luii i inosit y or mass (iMarconi fc Hunt 



l2003nillArg fc Rixl [2001' ' 
dark matter halo 



mass 



Ferrarese fc Ford I200E 
(jFerrarese fc FordI I200E 



early type stellar systems, one expects Mbh ~ 300 — 
10** M0 - although these correlations may c hange or 
break down at low mass (iGreene. Ho. fc Barth 2008) . 
This mass range is estimated by using the mea- 
sured Ursa Minor velocity dispersion (| Wilkinson et al.l 
120041 ) and by consid ering stellar mass-to-light ra- 
tios of ~ 1 — 3 (jZaritskv. Gonzal ez, fc Zabludofil 
20061) and a dark matter virial mass ~ 10^ Mo). 



Maccarone. Fender, fc TzioumisI ()2005l ) estimate an ex- 



http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg. pi 
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Fig. 1. — Low-energy (< 2.6 keV; above) and high energy (> 
2.6 keV; below) XISO+3 spectra. Shown arc total (red symbols), 
NXB (green symbols), and total— NXB (black symbols). The best- 
fit CXB-f GXB model to the NXB-subtracted spectra is shown by 
the black histogram. The best-fit was obtained by simultaneously 
fitting the 0.45-10.4 keV XISO+3 and 0.45-7.6 keV XISl spectra; 
the extension of the model, as well as that of the data, is plotted 
over the full 0.3-12 keV bandpass. 

pected X-ray luminosity of 7 x 10'^'* erg s^^ for Mbh = 
2.3 X 10"* Mq. We plot the Suzaku spectrum for an in- 
dex 1.7 power law with Lx = 3 x 10^'' erg s~^, and 
an optimistic estimate of the unresolved stellar emission, 
along with the background-subtracted XISO-t-3 spectrum 
of Ursa Minor in Figure 3. 

Typically, in Suzaku extended source sp ectral analysis 
(iMatsush ita et al.l l2007t iT akei et ani2007t iTawara et all 
[2008; Tok oi et all 12008: Sa to et al.ll2008D . the estimated 
NXB is directly subtracted and the source model often 
extended to include the GXB and CXB (perhaps with a 
subset of GXB/CXB parameters fixed based on previous 
observations). Offset or "blank-sky" fields may be em- 
ployed to directly subtract or fix the parameters of these 
components. The NXB-subtracted spectrum is binned 
and best fits and confidence levels are derived from 
statistics. Many previous limits on sterile neutrino radia- 
tive decay line fluxes are derived employing variations on 
these procedures, adding an "extra" line component and 
computing the line flux that results in an increase in 
by some amount corresponding to the desired confidence 
level. The appropriate Ax^ may be chosen under the 
assumption of Gaussian statistics (in which case Ax^ is 
distributed as x^, where q is the number of parameters of 
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Energy (keV) 

Fig. 2. — Same as Figure 1 for the XISl spectrum. 



XIS03 Background Subtracted Spectrum 




-t).5 1 1.5 2 2.5 



Energy (keV) 

Fig. 3. — Background-subtracted XlSO-f 3 spectrum of Ursa Mi- 
nor (crosses) compared with optimistic estimates of intrinsic X-ray 
emission from unresolved stellar sources (black histogram) and an 
AGN associate with an intermediate mass (Mbh = 10* Mq) black 
hole (red histogram.) 

interest) or base d on simulations l|Abazaiian "eTaLllMl . 
iProtassov et al.l ([2002) and others have identified short- 
comings of this general approach, although some may be 
obviated by carrying out simulations. 

We developed and applied a carefully considered al- 
ternative approach that we believe is statistically well 
founded. We do not subtract the NXB and use the un- 
binned data that is more appropriate for searches for 
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TABLE 1 
Exposures Times and 
Spectral Counts 



Detector 


Time 


Counts 


XISO+3 


138900 


14170 


XISl 


69450 


10610 



a narrow, weak spectral line. Spectral fits must thus 
include an NXB component. Best fits are determined 
through minimization of the maximum likelihood Cash 
statistic (C-statistic; Cash 1979), as appropriate for un- 
binned and unsubtracted data that obeys Poisson statis- 
tics, using version 12.4 of the XSPEC spectral fitting 
package^ ^. 

A baseline model is constructed as follows. (1) We fit 
the generated NXB spectra (Section 2.2.2) over the full 
0.4-12 keV band with a model consisting of two power- 
law components plus 10 narrow Gaussian lines with en- 
ergies fixed at values given in the Suzaku technical de- 
scription^^. XISO-l-3 and XISl spectra are fitted inde- 
pendently at this stage. (2) Next, we partition the band- 
pass into five intervals (0.45-0.65 keV, 0.65-1.3 keV, 1.3- 
2.6 keV, 2.6-5.2 keV, 5.2-E^a:r keV, where Y^max =10.4 
keV for XIS0-H3 and 7.6 keV for XISl) and refit each 
NXB spectral segment with the model described above 
(including only those lines that fall within the segment 
under consideration), allowing line energies and normal- 
izations to freely vary. If fits are improved by the inclu- 
sion of additional lines, by the replacement of Gaussian 
with Lorentzian line profiles, or by introducing finite line 
widths, these adjustments are implemented. The parti- 
tioning results in two line-free (0.65-1.3 keV, 2.6-5.2 keV) 
and two line-rich segments (1.3-2.6 keV, 5.2-Ejnax keV) 
with total counts fairly evenly apportioned among these 
four (see Table 1 for total bandpass counts for each de- 
tector). (3) All 10 XISO-l-3 and XISl NXB sub-spectra 
are then simultaneously fitted with line energies fixed 
at their locally determined values and one of the two 
power laws constrained to be the same for all spectra. 
This provides a phenomenological baseline NXB model. 
A baseline GXB+CXB model consists of an unabsorbed 
thermal plasma (apec) model, and variable abundance 
plasma (vapec) and power-law (plaw) models absorbed 
by the Galactic column density in the direction of Ursa 
Minor, as described in the previous section. The base- 
line NXB and GXB-I-CXB models are concatenated to 
form the total baseline model that we use as a starting 
point for our ultimate best fits. The best-fit vapec com- 
ponent require a sub-solar oxygen abundance - all other 
elements are consistent with solar abundances and fixed 
as such. 

The total Ursa Minor spectra are likewise partitioned 
into 10 segments (five each for XISO-l-3 and XISl). The 
baseline model is applied and then new simultaneous 
fits to the 10 spectra are found by minimizing the C- 
statistic. All NXB line energies and one NXB power-law 
component spectral index are now fixed; all constrain- 
able parameters for the GXB (two temperatures, one 

http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/prop_tools 
/suzaku_td/suzaku_td.html 



oxygen abundance, two normalizations) and GXB (one 
power-law index, one normalization) components may 
freely vary but are tied for all 10 spectral segments since 
they are best determined globally. GXB+CXB model 
components differ only by a constant multiplicative scale 
factor for XISl relative to XISO-l-3, to allow for resid- 
ual flux cross-calibration offsets. All NXB line fluxes, 
one power-law index and two power-law norms are al- 
lowed to independently vary in each spectral segment 
for each detector (allowing both power-law indices to 
vary independently did not improve the fits). Backfit- 
ting the best-fit GXB-I-GXB model components to the 
NXB-subtracted spectra described above yields accept- 
able fits; however, this is not the case when fitting the 
estimated NXB spectra with the best-fit NXB compo- 
nents. Evidently there are non-negligible (for our pur- 
poses) differences between the estimated and true NXB. 
These are apparent as residuals above 4 keV in the NXB- 
subtracted spectra (these residuals do not ruin global 
fits because of the relatively small number of > 4 keV 
source counts), and introduce systematic effects in the 
standard approach that we automatically account for 
in our analysis. When considering the inclusion of an 
additional (sterile neutrino) emission line, allowing the 
full range of NXB+GXB-I-CXB parameters of interest to 
vary provides a relatively conservative limit, while also 
utilizing unbinned, unsubtracted spectra enables us to 
consider energies for the extra line up to ^ 10 keV (in 
XISO-l-3). The issue of placing up per confidence limits 
on a quantity bounded from below (' Feldman fc Cousind 
[1998; Protassov et al. 2002) is mitigated since we can al- 
low line fluxes to go negative without the total spectrum 
falling below zero.^'^ 

2.2.3. Spectral Fits 

The best-fit NXB+GXB+CXB model to the Ursa Mi- 
nor spectrum constitutes the null hypothesis that the X- 
ray spectrum consists solely of background components, 
to which we wish to test the hypothesis that an addi- 
tional narrow"'^'' emission line - presumably from sterile 
neutrino radiative decay - may be present. We empha- 
size that the spectra are minimally processed - screened 
of spurious events, but extended to the maximum useful 
bandpass, and unbinned and unsubtracted. 

The final model has 62 parameters as follows: the 8 
GXB+CXB parameters described above (including the 
XISl/XISO+3 muhiplicative offset) and 54 NXB param- 
eters - one power-law index and two power-law normal- 
izations for each of the 10 spectra, nine emission line nor- 
malizations for each detector, plus two additional lines 
for XISO+3 with energies above the XISl maximum en- 
ergy of 7.6 keV, and two line widths for each detector. 
The best fit (Figure 4) has a minimized C-statistic of 
2863 for 2879 bins (2817 degrees of freedom); best fit pa- 
rameters are displayed in Tables 2 and 3 for GXB+CXB 

The "shce method", wherein the total spectrum is parti- 
tioned into slices corresponding to the appropriate energy resolu- 
tion, and the total counts in each slice us ed to derive an upper limit 
llRiemer-S0rensen. Hansen. &: PedersenI [20061) may provide a yet 
more conservative limit. This method, unnecessarily in our view, 
discards well-established characteristics of known X-ray emission 
components. 

Dwarf spheroidal dark matter velocity dispersions, 10 — 30 km 
s~^, are much less than the XIS energy resolution. 
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TABLE 2 
Best Fit CXB+GXB Parameters 



Energy (keV) 

Fig. 4. — Low-energy (< 2.6 keV; above) and high-energy (> 2.6 
keV; below) XISO-l-3 and XISl spectra; the XISl spectrum has the 
higher count rate. Each spectrum is divided into five segments, 
dehneated here by different colors, and all ten segments are simul- 
taneously fitted as described in the text. Histograms trace the 
best-fit model. 

and NXB parameters, respectively (in the former the 
XISl/XISO-l-3 factor was found to be » 1.1). Parame 
ters are consistent with previous work ([Matsushita et al 



2007t iTakei et alj l2007t ISato et al.l [20081 : IWerner et al 
2008f ) within their established variation on the sky as 
well as with the best fits to NXB-subtracted spectra of 
Ursa Minor described above. 

2.3. Sterile Neutrino Emission Line Flux Upper Limits 

The null hypothesis asserts that the Ursa Minor X- 
ray spectrum consists solely of background (NXB, GXB, 
CXB) components, with estimated parameters corre- 
sponding to the minimum value of the C-statistic (Cnuii)- 
We may add a narrow (sterile neutrino radiative decay) 
line and independently refit the augmented model to ob- 
tain Ciine(-Eiine) for a range of line energies i?iine- Given 
the large number of counts in our unsubtracted spec- 
tra, AC = Cnuii — Ciino(£^iinc) IS ideally expected to be 
distributed as the statistic with one degree of free- 
dom, and the con fidence level that a lin e is detected 
by P{xl < AC) (jKurczvnski et al.|[20M) . In the ab- 
sence of a significant detection, an a-upper limit line 
fiux is that which corr esponds to a value of AC such 
that P{xl < AC) = a (|Yaaoobl ll998l. 

We perform Monte Carlo spectral simulations in or- 
der to test how accurately these expectations are re- 



Model 


kT or r 


O Abundance 


Norm 


apec 


0.075 




0.0058 


vapec 


0.22 


0.21 


0.0016 


plaw 


1.45 




0.00087 



Note. — Temperatures in keV; 

oxygen abundance s relative to solar 

l lAnders fc Grevessd 11989 ): plasma model 
norms in units of J DVuenf^/ li-Kd? , 

where d is the source distance in cm and rie 
and nj, are the electron and hydrogen densities 
in cm"''; power-law norm in units of photons 
keV~^ cm~^ s~^ at 1 keV (corresponding to 
6.08 X 10~^ erg cm'^ 8"^ over the 2-10 kcv 
band). Norms correspond to a 20' circular re- 
gion, that is, 1256.6 sq-arcminute=0. 0001063 
str. 

alized. The parameters from the best-fit null hypoth- 
esis model to the observed data are used to generate 
fake spectra that are re-fitted with the null hypothesis 
model. These second generation fits a re used as seeds for 
subsequent Monte Carlo realizations (iButler et al.l 120051 : 
iMarkowit z et al.ll2006l:lPorauet et al.l l200/) that are gen- 
erated and fitted in a manner that matches the actual 
spectra and their best fits (same exposure times and in- 
strument responses, same assignment of fixed, tied, and 
free parameters). Fitting these simulated spectra with 
the null hypothesis model determines the minimized fit 
statistic CnuU- Each simulated data set is additionally 
re-fitted with models that include an additional line at 
an (non-fixed) energy E'line that is stepped through the 
analysis bandpass (at 25 eV intervals) to determine the 
minimum fit statistic Ciinc(£'iine) for each Sunc- 1000 
simulated spectra are generated for each hypothesis (null 
or null- plus- line- at-_Eii,ic for each inline)- The cumu- 
lative probability distribution P{< AC), where again 
AC = Cnuii — Ciine(-Eiine), yields an experimental esti- 
mate of the corresponding probability of rejecting the 
null hypothesis. The upper limit may then be esti- 
mated as described above, using the experimental, in- 
stead of Xi distribution. We find that P{< AC) matches 
P{< x\) only in spectral intervals with high total count 
rates and low NXB (i.e., in NXB-line-free regions). The 
joint distribution over all energies is well-approximated 
by P(< x^), and at no energies does P{< AC) fall below 
P{<xl) (Figure 5). 

We adopt AC = 10 upper limits. This corresponds to 
P{< AC) ~ 0.99 in the joint distribution; at no energy 
was P(AC < 10) less than 0.984 (or greater than 0.996). 
The AC = 10 upper and lower limits are plotted in Fig- 
ure 6. Since no lower limit is > 0, we infer that there is 
no line detected at > 99% confidence. To derive limits in 
the vicinities of NXB lines we fix their line fiuxes (oth- 
erwise, in fits and simulations, these are free to vary) at 
their best-fit values. The resulting limits are very nearly 
equal to the AC = 10 uncertainty in the NXB line fiuxes 
and, as they neglect any systematic errors in the NXB 
line strengths, will underestimate the upper limit of any 
additional emission line with precisely the same energy 
as the background features. The limits are tightest in 
spectral intervals where lines from the NXB and GXB 
are sparse or absent (especially '^2 — 5 keV). When we 
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Fig. 5. — Cumulative probability of AC = Cnuii — Cunc from 
simulations. The solid curve shows the sum over all energies, while 
the broken curves shows the full range spanned by simulations at 
individual energies. Crosses, boxes, and triangles show the 68, 90, 
and 95% cumulative probabilities for the Xii xi' ^''^'^ xi distribu- 
tions, respectively. 



-3.5 




2 4 6 

line energy (keV) 



Fig. 6. — Log of AC = 10 line flux upper limits in photons cm~^ 
s~^ (solid histogram). Also shown is the log of the absolute value 
of the lower limits (that are all negative; broken histogram). 

compare these limits to those more conventionally de- 
rived from the binned, background-subtracted spectra, 
we find that the AC = 10 limits roughly correspond to 
Ax^ = 7 in the < 4 keV region where the NXB is small 
or cleanly subtracts out. 

3. LIMITS ON STERILE NEUTRINO PARAMETERS 

The integrated line flux from sterile neutrinos decaying 
at rate Fgt = 10~^^r_27 s""'^ from a dark matter mass in 
projection of lO^Mj M© at a distance lOOdioo kpc is 

^line = 7.5 X 10-15 r_27/stA/7droo ^rg cm-^ s-\ (2) 

where f^t is the dark matter fraction in sterile neutrinos. 
The corresponding number flux of E-^ — mst,/2 photons 
from decays of particles with mass rrigt that make up the 



(3) 



!roo (4) 



Fig. 7. — Solid (broken) curve shows the 90% c onfidence lim- 
its of the least (most) compact mass profile from IStrigari et al.l 
(2008a ). Al so shown are estimates, at various ra dii from 
Gilmore et al,| ||2007) (G07), Wu (2 0071) (W07'). iStrigari et al.l 
(2007b) (S07b), Penarrubia, McCoi machie. fc Navarrol (I2008IT 
(P08), and .Strigari et al, ( 2008b ) (S08b). 

total mass 10''M7/st Mq is given by 

^ii„e = 9.4 X 10-^ r_27 /st A/y^roo f— ) 

cm"'^ s""'^, where we have used the expression for T^t 
introduced in equation (1). 

One can use this relation to set limits on sterile 
neutrino mass and mixing angle using the photon 
flux limits shown in Figure [H The dark matter halo 
mass distribution in Ursa Minor may be estimated 
using stellar dynami cs. Strig ari and collaborators 
(IStrigari et all l2007allbl . l2008allbl ) applied maximum 
likelihood analysis to constrain the parameters of 
a generalized Nava rro- Frenk- White (NFW) profile 
(jPiemand et al.ir2005f ). Based on their confidence region 
for an NFW model (inner d ark matter density slope of 
— 1) in [Strigari et al.' (2008a), and the general integrated 
constr aints (Gilmore et al. 2007; Wu 2007; Strigari et ap 
2007U fP cnarrub ia. McConnachie. fc Navarrol |20^ 
Strigari et al.,.2008bl ). the projected mass of Ursa Minor 
within 20' (~ 0.4 kpc at our adopted distance of 69 kpc; 
Grebel, Gallagher, & Harbeck 2003) is well determined 
from t he projection o f the N FW model (Bartclman^ 
199fit iGolse fc KneibI [2001 : lYang. Mo. fc Kauffma^ 
2003f) to be ~ 6 X lO'^ Mq. This is consistent with 
the weighted average derived from the estimates in the 
papers quoted above (Figure 7; these adopt various 
assumptions about the cuspiness of the profile) that 
has a variance of ^ 25%, althoug h the ra nge of models 
formally allowed in [Strigari et all (|2008a[ ) spans factors 
of ^ 2 in either direction at that radius. The correspond- 
ing excluded dark matter decay rate Fst (for fst = 1), 
derived from the AC = 10 line fiux upper limits, is 
shown in Figure 8 for line energies 0.45 < ekoV < 10.4. 

If one assumes that all the relic sterile neutrinos are 
produced via the DW mechanism, the absence of line 
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emission from Ursa Minor provides an upper limit on the 
mixing angle as a function of mass that is independent 
of whether there are additional production mechanisms 
(which could only increase the abundance). The flux 
from equation (jj]), calculated for a given value of fst, 
should not exceed the photon flux shown in Figure [6l In 
a given cosmological scenario, if the sterile neutrino relic 
abundance is a monotonic function of mass and mixing 
angle, fst = p{sm^ 9, rUst), then the region allowed in the 
sin^ 9 — plane is given by 

sin^ 9max,p = min(sin^ 9,naxs, sin"^ 9max,2), (5) 

where sin^ 9 max. i is determined from fst — 
p(sin^ 0, TOst) < 1 and sin^ 9max,2 from the require- 
ment that the flux in equation ((4|) not exceed the 
observed flux shown in Figure 6. 

We derive two different limits that provide answers to 
the following questions: (1) whether the existence of a 
sterile neutrino with a given mass and mixing angle is 
consistent with standard cosmological history and (2) 
whether sterile neutrinos of a given mass and mixing an- 
gle can account for 100% of dark matter. These two 
constraints are, of course, different. 

For the first limit, we take fst — PDw(™st,6') to be 
a function of mass and mixing angle and we assume 
the DW production mechanism, which gives the minimal 
abunda nce of sterile neutrinos in standa r d Big Bang cos- 
mology. lAsaka. Laine. fc Shaposhnikovl (|2007h provided 
fitting formulae for sterile neutrino abundance produced 
by the DW mechanism of the form 

PDw(sm 9, mst) = a j j , (6) 

as well as pairs (a, (3) for average, minimum, and maxi- 
mum production that account for hadronic uncertainties. 
The excluded region in the mst — 9 plane, for production 
dominated by the DW mechanism, is given by equation 
(4), with 

fsm^9max.i\ -1 / rUst \-f^ 

\ 10-1° j^" ' 

and 

Sin^ 9ynax,2 = (siu^ 9ynax,l X sin^ 9,naxfl[^]Y^'^, (8) 

where sin^ ^maa, o[l] is determined by setting fst = 1- 
The resulting exclusion region is shown in Figure [9] for 
the mini mal rate of DW production cons i stent w ith the 
results of lAsaka. Laine. fc S haposhnikovl (|2007D . which 
provides the most conservative constraints. This region 
is excluded regardless of the physics responsible for mass 
generation of neutrinos or any other physics beyond di- 
rect mixing between sterile and active neutrinos. 

For the second kind of limit, we determine the part of 
parameter space for which sterile neutrinos can account 
for all the cosmological dark matter, while still being 
consistent with Suzaku observations. Here we set fst = 1 
without reference to any specific production mechanism. 
The excluded region (for fst = 1) corresponding to the 
AC = 10 upper limits on Fst/st shown in Figure 8 is 
delineated by the solid line in Figure [9l Sterile neutrinos 
occupying the parameter space to the right of the solid 
line cannot make up all of the dark matter (although 
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Fig. 8. — Excluded dark matter radiative decay rates (hatched 
region) as a function of energy averaged over the inner 0.4 kpc of 
Ursa Minor, assuming that 100% of the dark matter is composed 
of sterile neutrinos (/st = 1). 
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Fig. 9. — Sterile neutrino parameter space to the right of the 
solid curve is excluded by the Suzaku observation of Ursa Minor 
if dark matter is solely composed of sterile neutrinos produced 
by some (unspecified) mechanism. The solid exclusion region is 
model-independent, based only on the assumption of the standard 
cosmological history below the temperature of a few hundred MeV, 
when the DW production by neutrino oscillations takes place. 

such a particle may exist in nature, in contrast with the 
previous limit). 

One can ask which dark-matter particle mass in the 
form of sterile neutrinos, for fst — 1, produced solely 
by the DW mechanism is consistent with Suzaku obser- 
vations. Here one can set an upper limit on the sterile 
neutrino mass: rUst < 2.5 keV. However, the DW mecha- 
nism with /st = 0.1 is not ruled out at any energy covered 
by the Suzaku spectra ('^1 — 20 keV). This is shown in 
Figure 10, where our line flux upper limits are compared 
with the DW predictions. 

3.1. Comparison to Previous Work 

Because different criteria for estimating line flux up- 
per limits are used, a precise comparison with pre- 
vious results is difficult. Our limits are about a 
factor of 3 better than those fo r Ursa Minor in 
iBovarskv. Nevalainen. fc Ruchavskivl ()2007[ ) over the 
overlapping bandpass. The constraints shown in Fig- 
ure 9 are compara ble to the combined constraints from 
the Milky Way (|Bovarskv. Nevalainen. fc Ruchavskivl 
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Fig. 10. — Line flux upper limits from Figure 6 re-plotted 
with the /st = 1 (upper curves) and fst = 0.1 (lower curves) 
DW predictions for maximum, ave rage, and minimum production 
llAsaka. Laine. &c Shaposhnikovll2067l '). 



I2007f) and M31 (|Bovarskv et all |2008[ ). and slightly 
tighter than those from the unresolved Chandra CXB 
([Abazaiian et al.1 [2007') an d from the "bullet" clus- 
ter, IE 0657-56 CB oyarsk yTRuchavskiy. fc MarkevitchI 
|2008() . lAsaka. Laine. fc Sha poshnikovl (|2007l ) compiled 
and plotted various X-ray limits on the combination of 
TTist and 6. We attain a comparable level of sensitivity 
over the entire 1-20 keV mass range with a Suzaku ob- 
servation of a single object. 

4. DISCUSSION, CONCLUSIONS, AND FUTURE 
PROSPECTS 

With a single Suzaku observation of the Ursa Minor 
dwarf spheroidal, we have derived limits on the mass 
and mixing angle of the sterile neutrino dark matter 
candidate that are comparable to or better than previ- 
ous constraints. These limits are illustrated in Figure [9] 
for the case where sterile neutrinos comprise all of the 
dark matter regardless of how they are produced in the 
early universe and, in the more general case, when their 
abundance is not assumed, but is calculated under the 
assumptions that minimize it. The best estimate of the 
total projected dark matter mass within 400 pc of 6 x 10^ 
Mq is adopted here. If solely produced by the DW mech- 
anism, sterile neutrinos cannot constitute 100% of the 
dark matter in Ursa Minor unless they are less massive 
than 2.5 keV. 

Because of its orbit, Suzaku has a lower and more stable 
background than either the Chandra or XMM-Newton 
Observatories, and its spectrometers have sharper energy 
resolution. Our approach in deriving limits utilizes the 
C-statistic to evaluate a very general model fit to the data 
with minimal pre-processing and no background subtrac- 
tion. Thus, there arc no a priori assumptions about the 



accuracy of the NXB estimate that we only use to set line 
energies and shapes, or whether an extra line feature lies 
above or below the total continuum. Our final limits 
are those that correspond to a change in the C-statistic 
that occurred in 1% of our Monte Carlo simulations. 
These factors may account for why our constraints rep- 
resent only modest improvements over previously pub- 
lished limits despite the more optimal characteristics of 
Suzaku. As the NXB characterization improves, it may 
be possible to tighten upper limits in the rrist > 2 keV 
regime; at lower masses, the ubiquity of the CXB will re- 
main a barrier until significantly higher energy resolution 
is attained - e.g., with Astro-H^^ . 

Since there is some uncertainty as to the dynamical 
state of dwarf spheroidals on an individual basis. X-ray 
constraints obtained from additional systems are war- 
ranted (e.g., our paper - in preparation - on the Suzaku 
observation of the Draco system). The Sloan Digital Sky 
Survey (SDSS) has approximately doubled the number 
of galaxies identified with the local group population and 
revea led a new sub-populati on of faint Milky Way satel- 
lites ([Belokurov et al.l l 2007D. including some with very 
high mass-to-light ratio (jStrigari et all l2008al |bl) that - 
with the right observational setup - may provide the best 
targets for investigations of this type. 

Finally, predictions of the Lya forest power spectrum 
in fully self-consistent large-scale structure simulations 
with sterile neutrinos can provide complementary con- 
straints in the form of a lower limit on the particle mass. 

Sterile neutrinos remain a viable dark matter candi- 
date. For example, a sterile neutrino with a mass of 
mst ~ 4 keV and sin^ 6* ~ 4 x 10"^" could account 
for all of the dark matter (or ~ 15%— 75% if produced 
solely by the DW mechanism) and exp lain pulsar kicks 
(|Kusenko. Mandal. fc Mukherie^ 120081 ) without violat- 
ing the line flux upper limits. 
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TABLE 3 

NXB Line Energies, Best-fit Fluxes, and Flux Upper Limits 
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Note. — id entification s, when available, from Table 7.2 of "The Suz aku Technical Description" 

ihttp: / /heasarc.gsfc.nasa.gov/docs/suzaku/prop_tools/suzaku_td/suza ku_td.html| ; line energies in keV, lines 
luxes in iO""^ photon cm^'^ s~^, upper limits correspond to AC — 10. 
^ Lorentzian^ finite width 



